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ABSTRACT 



Experimental investigations were undertaken to determine the collapse 
pressures of relatively thick circular cylindrical shells subjected to 
hydrostatic pressure. A total of twenty-three fiber glass reinforced 
plastic cylinders with thickness to diameter ratios ranging from 0.C5 to 
0,09 was tested to determine collapse pressures and the effective modulus 
of elasticity of the material. The length to diameter ratio was kept con- 
stant. Results were plotted non-dimensionally with P/E as the ordinate 
sind h/D as the abscissa. The results were compared with available 
theoretical instability formulas and were in very good agreement with 
theory for ratios of h/D up to about 0.07. Above this value of h/D, 
experimental collapse pressures were lower than those predicted by 
theory. 

It was concluded that the von Mises equation [equation (6), Ref. 22] 
is the best instability equation for cylinders with ratios of h/D greater 
than 0.05. It was also concluded that for almost all practical cases, 
the thin shell theory is sufficiently accurate for predicting instability 
failures. 

Derivations for the following are included: (l)Expressions for 
calculation of the moduli of elasticity in the circumferential and 
longitudinal directions; (2)An expression for the shear modulus and 
for the effective modulus of elasticity of an orthotropic material; 

( 3 ) A determinant expressing the critical pressure of simply supported, 
orthotropic, circular cylinders subjected to hydrostatic pressure; and 

(4) An expression for calculating the bending stress in a simply supported 
circular cylinder subjected to hydrostatic pressure. 
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NOTATION 



D 

D 

D 1 

D2 

E 

Ee 

Ex 



Flexural rigidity, D - 



cA’ 

iz(/-v l ) 



«* 



Mean diameter 
Inside diameter 
Outside diameter 
Modulus of elasticity 

Effective modulus of elasticity for anisotropic material 
Modulus of elasticity in the longitudinal direction 
Ee, Modulus of elasticity in the circumferential direction 



G 

h 

K 1 

k 2 

L 

^oa 

m 

m* 

Mx 

M# 

M x ^ 

M*x 



Shear modulus, G = tt-, -? 

* 2 ( l+v) 

Thickness of shell 



Non-dimensional term from thick cylinder stress formula, 

zRl 

Ki = rF W 

Non-dimensional terra from thick cylinder stress formula, 

v. _ 

K2 'W^' K ‘ 

Effective length 
Actual length 



Number of waves or lobes along effective length of cylinder at 
time of collapse 

Strain-pressure ratio, circumferential strain gage 
Strain-pressure ratio, longitudinal strain gage. 

Axial bending moment per unit of circumference 
Circumferential bending moment per unit of length 

Twisting moment per unit of length on element cut out by generators 

Twisting moment per unit of circumference on element cut out by 
plane perpendicular to axis 



«* 

In the notation used in this paper, the subscript "1" refers to the 
inside of the cylinder, and the subscript "2" refers to the outside of 
the cylinder. 
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n 

Nx 

Np 

N** 

N^x 

P 

Qp 

Q* 

r 

R 

R 1 

r 2 

u 

V 

w 

z 

ot 

r 

€x 

60 

0 
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Number of lobes or waves in a complete circumferential belt around 
cylinder at time of collapse 

Normal force per unit of circumference in axial direction 

Normal force per unit of length in circumferential direction 

Shear force per unit length on elements cut by generators 

Shear force per unit of circumference on elements cut by planes 
perpendicular to axis 

Pressure 

Radial shear in the xz plane 
Radial shear in the yz plane 
Radial coordinate 
Mean radius 
Inside radius 
Outside radius 
Axial displacement 
Circumferential displacement 
Radial displacement 

Perpendicular distance from mid-surface of shell, positive inward 

i 

Non-dimensional terra in buckling equations, oi = 3" pi 
Shearing strain 

Shearing strain along i-plane in j -direction 
Normal strain in axial direction 
Normal strain in circumferential direction 
Normal strain in radial direction 
Cylindrical (angular) coordinate 

Tr D 

Non-dimensional terra in buckling equations, \= m Y77 
Poisson's ratio 



v 



O'x Normal stress in social direction 
o> Normal stress in circumferential direction 
Shear stress along i-plane in j-direction 
4> Cylindrical (angular) coordinate 

'X, x Change of curvature of generator 

Change of curvature of circumference 
%xj< Twist of element due to M x p 
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INTRODUCTION 



The absence of experimental results and elastic instability analyses 
for thick cylinders has been pointed out by Wenk^. On the basis of this 
information and consultation with Mr. John Pulos of the Structural Mechanics 
Laboratory, David Taylor Model Basin, the authors decided to undertake an 
experimental investigation of the elastic instability of relatively thick 
cylinders subjected to external hydrostatic pressure. Thus, the primary 
purpose of this thesis evolved; which is the experimental evaluation of avail- 
able theoretical instability formulas as applied to cylinders with thickness 
to diameter ratios greater than those heretofore tested. From the results 
of the proposed experiments it was anticipated that the following question 
posed by Dr. Wenk (see "Discussion” of Ref. 1) could be answered: 

"At what ratio of thickness to diameter are errors in the thin-shell 
approximation unacceptable?" 

A prime requisite prior to undertaking the thesis was the making of a 
careful survey of the literature to determine the extent of theoretical and 
experimental work. The authors conducted 3uch a survey both prior to 
beginning the experiments and on a continuing basis throughout the period 
of the thesis work. The results of this search are presented in the section 
entitled "Literature Survey." 

It was also necessary to design and construct the experimental appa- 
ratus and the models. The materials and methods of design, together with 
the procedures employed, are indicated in the body of this paper. Since the 
models were fabricated from glass- reinforced plastic, the elastic formulas 
required in analyzing the data were modified to account for the difference 
in the values of the principal moduli of elasticity. Derivations for the 
expressions used are presented in the appendix. 

1 

References are listed beginning on page 67. 



The thesis begins with a section on concepts and assumptions employed 
in the theory of thin shells. This study was made in order to gain an 
insight into the differences between a thick and a thin shell theory. 
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CONCEPTS AND ASSUMPTIONS IN THE THEORY OF THTN SHFTT.q 



A thin shell may be defined as one for which the ratio of the thick- 
ness of the shell to the principal radii of curvature is small compared 
to unity. The problem of classifying a shell as thick or thin therefore 
reduces to the determination of a value of h/R below which solution errors 

p 

will be within a prescribed limit. Novoxhilov has stated that for or- 

technical calculatxons it is admissable to have a relative error 

of 5%, and that for thin shells this corresponds to a maximum h/R = 0.05. 

3 

Timoshenko (p.354) has suggested a value of h/R = 0.10 as the dividing 
point for thin curved bars. 

There are three principal assumptions employed in the theory of thin 
shells in addition to those made for most problems in elasticity. They are 

(1) The normal stresses acting on planes parallel to the middle 
surface may be neglected in comparison with the other stresses. 

(2) The normals to the undeformed middle surface are deformed with- 
out change of length into the normals to the deformed middle 
surface. 

(3) The thiokness h is very small compared to the least radii of 
curvature (i.e. h/R«l). 

Assumptions (1) and (2) essentially reduce the problem to the case of 

plane stress. Consider the generalized matrix form of Hooke's law (see 
4 

Borg , p. 59). 

T = Al^E^ + ^ /U/J 

where 

T = the stress tensor 

- the strain tensor 

1^= the trace of the strain tensor 
the unit, 3x3, matrix 

- constants 

E £ 

^ ( / + -%/)(/-£ -i/) y° ^ 
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Now at the inside surface of the cylinder (i.e. at z = h/2) the 
stress component <r z is zero. Then, if the shell is thin, this compo- 
nent can be assumed to be zero throughout the thickness of the shell 
since it is small compared to the stresses c r# and a ~ x . 

Assumption (2) is equivalent to assuming that the shearing strains 
Y xz and Yq z are zero, for the definition of shearing strain is the 
change in value of an originally right angle from the unstrained state. 
The assumption is that this right angle is preserved during deformation. 

Therefore, letting cr_ = V' =Y' be zero, the tensor equation 
becomes 



r 


&KQ 


<Xz \ 


' 


0 


i <^e 






= A(6^6 e + € x ) 0 


/ 


\<r xz 


&OZ 


of 


\o 


0 



0 + ZyU 




0 

0 



Expanding the equation gives 

= A (6* + + £*) + ZyU € k 

a~ e = + ZyCt £ 0 

<Jl — O - A (€< *■ + € z ) 

Ycq Y. 0 

ft*. ~ z — O 



Eliminating £ z from the first two equations by use of the third gives 
the desired stress-strain relations 



£ 

= ~Z (t -ts) 6 

Furthermore, assumptions (1) and (3) are utilized in reducing the 
strain components to a usable form (see Wang'’, pp. 335-340). When this 



- A- - 



^Xa)] 



is done, the equation for cr^ becomes 
E r 

/ — JS *- fixO + O _ Z C X * +■ 

where c> h inf. 

“7x '*•* = 

__ J- in’ _ nr v _ _l ^ ( nc +. jl i *£ ) 

^ 0O <3- cJ O <3- 0 ^ \ £«_ ^ & / 

The expressions for the other stress components are similar in degree 
of complexity. 

Finally, to obtain the relations between the stress resultants 
and couples, the following expressions must be integrated. 





Again, assumption (3) is used and z/R© is neglected. Then, N xe = N ex , 
Mjjq = Mgjj, and the expressions may be integrated without difficulty. 
For example, the expression for N x becomes 



N. 



Eh f 

/ - y * L ^ 



i-C 



■f y 






djyT 



^j- 



)] 



This equation and the remaining ones obtained in a similar manner are 
known as Love's first approximation. 

It is possible to derive expressions for the stress resultants and 

couples without resorting to the use of assumptions (1) - (3). Naghdi^, 

7 

using a variational method due to Reissner , has made such a derivation. 

In order to gain an appreciation for the simplicity of the above expression 
and of the complexities involved in a thick shell theory, the more exact 
expression for N x is given below. 



-S- 



N k =[tT^X Uo "'feo]} + *■ {-7^jk (x* -Xe) 

X 



v X* 



e ^ ' '<•* • ] 

Re J 



__^ L 

/ - >-" z. 



[_(**■)*"£ + - -Ji] j 



Presumably, these more exact relations could be introduced into the 
equilibrium equations, and, provided a solution could be effected, a 
thick shell analysis for the instability problem would become available. 
An analytical approach to a thick shell instability analysis has not 
been attempted in this thesis, however. 
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LITERATURE SURVEY 



The theoretical and experimental analysis of the buckling of 
cylindrical shells subjected to hydrostatic pressure has attracted the 
attention of numerous investigators dating from experiments conducted by 

O 

Fairbairn in 1858. All investigators employed the theory of thin shells 
in arriving at theoretical formulas and most experiments were carried out 
on models with relatively thin walls. Tables and descriptive information 
are presented in this section to indicate the extent of previous theoretical 
and experimental investigations. Particularly helpful in the making of 

q 

the literature survey were two reports by Nash which contain a listing 
of 2339 papers and books relevant to experimental and theoretical work 
on shells and shell-like structures from 1828 through 1956. 



Results of Previous Experimental Investigations 

Table 1 is a chronological listing of various investigators together 
with the range of geometries involved and other pertinent data. Since the 
first analytical solution to the instability problem was not available 
until 1913, all investigations prior to this time were made principally 
to determine empirical design formulas. Those investigations made in later 
years were carried out to test the validity of analytical solutions either 
in general or as applied to a particular problem. 

With the exception of Southwell, all experimental investigators prior 
to 1929 restricted the loading on the cylinder to a radial pressure only. 
The axial thrust was eliminated by the use of a relatively stiff rod which 
connected the end plugs. The use of radial pressure only resulted from the 



absence of analytical formulas based on combined loading, and from the 
fact that one of the principal applications of the results was to the 
study of boiler flues, Southwell did use a combined loading to simplify 
the apparatus. He neglected the effect of axial loading, however, in 
comparing test results with his analytical approach. 

Most investigators after 1929 employed models subjected to hydrostatic 
(combined radial and axial) loads. An analytical solution to the problem 
of the cylinder under hydrostatic pressure was available by this time, and 
the submarine presented an immediate application. 

From columns 4,5, and 6 of Table 1 the following is noted: Most of the 
models were relatively small in diameter, as was the h/D ratio. In those 
cases where the h/D ratio was greater than 0.05 the model was either so 
long that the values of the critical pressure approached that of an infinite 
tube, or they were in the range where yielding of the material occurred 
rather than an instability failure. 

The type of support used by the various investigators is indicated 
in column 7. In general, this information was procured from sketches of 
the models as prepared for testing. A clamped type of support is indicated 
only where the author concerned specifically labelled the test condition as 
such. It is again noteworthy that no attempt was made to clamp the ends 
until 1941^ when an analytical solution became available. Also, as indi- 
cated by Cook}^ simply supported ends can be approached experimentally 
much more closely than fixed ends. 

Column 8 indicates the maximum collapse pressures achieved by the 
various investigators; column 9 the number of models tested; and column 
10 the variety of materials which have been utilised in the construction 
of models. 

Several general observations can be made based on the survey of 
experimental results found in the literature: 



(1) Scale effect is not important provided reasonable care is taken in 
fabricating the models. After testing 100 models representing 
the strength hull of a submarine, Windenburg^ concluded that 
scale effect was not ^evident. Some the earlier experiments 

have not given good results when compared to theory because 
of the inadequacy of materials and unsound fabrication techniques. 
Most of the results obtained in the last two decades have given 
very good correlation. 

(2) Although the assumed end condition is taken as either simply 
supported or as clamped, neither condition can be precisely ob- 

11 

tained in practice. Based on tests of 258 brass models, Cornell ^ 
concluded that the end plugs inserted into the ends of the tubes 
offered a restraint that tended to reduce the effective length. 

He found that to get close agreement of his test results with 
theory, it was necessary to use an effective free length which 
differed from the actual free length by to (Because of 

the small value of such a correction, the authors assume the 
measured free length to be the same as the effective free length 
for the tests reported in this thesis.) 

(3) Several investigators have conducted carefully controlled experi- 
ments to determine the influence of other factors, such as eccen- 
tricity, on the collapse pressure. Sturm^ was probably the first 
to investigate the problem of out-of-roundness from an experimental 
and analytical point of view. In recent years, Cleaver^ has 
conducted very thorough, extensive, and carefully controlled 
tests on 530 models subjected to radial pressure. He lists the 
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following factors as affecting the collapse pressure: 

a. Eccentricity of bore relative to the external surface. 

b. Variations from true circular shape (out-of-roundnesa) . 

c. Ovality of the specimen (the difference in the maximum and 
minimum diameters divided by the nominal diameter) . 

d. Variation of material properties 0 

From the results of his tests Sturm concludes that eccentricity of bore 
is the predominant manufacturing vuriable affecting collapse strength, 
but for eccentricities within the limits imposed by current specifications, 
its effects are small and can be neglected for practical purposes. 

Variations in circularity exert no measurable systematic influence, and 
provided they are within current specifications (i.e. for commercial tubes), 
they too may be neglected. 

The authors point out in passing that for a well-made, simply 
supported model the above effects tend to cancel the effect of any 
end restraint that might be present. Therefore, neglecting these effects 
may be justified when attempting to correlate experimental and theoretical 
results. 



Results of Analytical Investigations 

Table 2 presents a summary of results of principal theoretical 
investigations. The method of approach to the problem is indicated. 

The differences in the various formulas can be accounted for by 
the inclusion or omission of certain higher order terms in the differ- 
ential equations of equilibrium, by the method of approach (i.e, equil- 
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ibrium, energy, Icinetic ), and by the degree of simplification of the 
final result. 

All of the more exact formulas are found to have the form 

p/E = (h/D) 3 + F 2 (h/D) 

where F, and F_ are functions of 
1 2 

(a) The length to diameter ratio. 

(b) Poisson's ratio. 

(c) The number of lobes in a complete cicumferential belt at the 
time of collapse. 

(d) The number of lobes along a generator of the cylinder at the 
time of collapse. Von Mises^* has shown this quantity to be 
unity for minimum buckling pressure. 

In order to facilitate comparison of the buckling equations with 
test results, values of F^ and F 2 for the model geometry and material 
under consideration in this paper have been calculated for the various 
buckling formulas. These values are given in Table 3. Curves for p/E 
versus h/D are also included for these cases. (See Fig. 13). 

In addition to the formulas for the buckling pressure of cylinders 
under hydrostatic load, it is of interest to examine formulas for buck- 
ling pressure under uniform radial load only. Table 4 gives some of the 
principal radial load buckling equations. A comparison of the radial 
load formulas with the hydrostatic load counterparts indicates that the 
following relationship is approximately true 

2 

buckling pressure under hydrostatic load = n - 1 

. . ■ ■ ... . - . ■ M . g 

buckling pressure under radial load vr - 1 + A 

2 

Thus as A approaches *ero (L/D « ) , the formulas for hydrostatic 
load should yield the same results as the formulas for radial load. 

For the case where A = 0, n becomes two, and all of the more exact 
hydrostatic and radial load buckling pressure equations reduce to the 
Bresse equation (Table 4) for a tube of infinite length. 
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TEST APPARATUS. MODELS. AND TEST PROCEDURE 



Test Apparatus 

The test apparatus consists of a pressure vessel, pump, tubing, 
pressure gages, and equipment required for the taking of strain gage data. 
The arrangement and details of the apparatus are shown in Figures 1 through 
5. Specifications which were drawn up for the test apparatus are given 
below. 

Size of chamber: 5" I.D., 13" depth 

Working pressure: 10,000 psi 

Penetrations: One fitting for pressure tubing, one fitting for 

four wires from strain gages on the external 
surface of the model, one fitting for l/4" tube 
to vent model to the atmosphere and to lead out 
internal strain gage wires. 

Pump: Hand operated pump with self-contained check valve and 

reservoir— the type commonly used with hydraulic jacks. 

Pressure gages: One 5000 psi 8" Bourdon gage graduated in in- 
crements of 50 psi, and one 15,000 psi 6" Bourdon 
gage, graduated in increments of 100 psi. 

Strain Indicator; Baldwin, Type N 

The pressure vessel was constructed in accordance with a basic design 
drawn up by Mr. Kenneth Horn of the David Taylor Model Basin. The design 
was modified slightly by th« authors, to suit the purpose required. The 
following calculations show the factors of safety, for a working pressure 
of 10,000 psi. The material is low carbon steel, the yield strength of 
which was assumed to be 30,000 psi. Formulas used are given in reference 
(39). 



-/?- 



t = wall thickness, inches 



,25 for walls i 
,125 for bottom) 

d = inside diameter of chamber, inches, = 5.00 



= f 2 * : 
* (2.: 



P = chamber pressure, psi 

S = maximum stress, psi, = 30,000 

R = inside radius of chamber, inches, = 2.50 

F.S. = factor of safety based on yield strength of 30,000 psi 
Based on circumferential stress in wall: 



P = = (30,000) (2.25) = 17,532 psi 

R + 0.6t 2.50 + (0.6) (2.25) 

F.S. = 1L122 = 1.75 

10,000 

Based on longitudinal stress in wall: 

p = -2§1— = Liiuo^ogou^ = 84 375 i 
R - 0.4t 2.50 - (0.4) (2.25) * P 

p 5 = ^4 i375 ~ 8.44 

* 10,000 

Based on stress in bottom: 

p - . st2 , ** = ( ?0|00g K ?-l^) 2 = 33,449 p ,i 

P d* 2 (0.162) (5.00)^(0.162) P 

ps - 22 1 44? = 3.34 

* 10,000 * A 



* Ref. (39), para. UG-27(c)(l) 

»* Ref. (39), Section VIII, para. UG-34(1). 

Para. UG-36(3) states: 

"Single openings in vessels not subjeot to rapid fluctuations in 
pressure do not require reinforcement other than that inherent in the 
construction under the following conditions: 

....(b) Threaded, studded, or expanded connections in which the hole 
cut in the shell or head is not greater than two-inch pipe size." 
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The pressure vessel was proof tested to 11,000 psi on 30 August, 
I960. It was raised to this pressure three times, and strain data were 
taken on the outside surface using a circumferential strain gage. This 
test was made prior to the making of the penetrations for fittings "A" 
and "B" (Fig. 4). The strain data taken plotted linearly from zero to 
11,000 psi. 

This apparatus should not be used with a working pressure in excess 
of 10,000 psi, unless the pressure vessel is again proof tested at the 
higher pressure with strain gage measurements taken in regions of high 
stress. In any case, it is recommended that a pressure of 15,000 psi 
not be exceeded, and that for working pressures in excess of 10,000 psi 
the pressure vessel be placed in a pit or behind a suitable barrier. 

All tubing and fittings used should be rated at the maximum working 
pressure employed. 

Persons making further use of this apparatus are warned that 
valve ”G" . Fig . ^ , i£ rated at 3000 psi. It is recommended that this 
valve be replaced by one rated at 10,000 psi or higher. As a safety 
precaution in performing the tests for this thesis, a 1/8” steel plate 
was attached to the table between the operators ^nd fittings "D" through 
"H" (Fig. 4). The plate can be seen in the photograph, Fig. 5. 
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FIGURE 5 

Photograph of Assembled Apparatus 
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Model Dealgo 

Design of the models was based on the following considerations; 

(1) Maximum working pressure of test apparatus; 10,000 psi 

(2) Dimensions of test chamber; 5" I.D. and 13" depth 

Since the purpose of the thesis is to evaluate the instability formulas 
for relatively thick-walled cylinders, it was necessary to design the models 
to fail by instability, and to have a relatively large thickness-diameter 
ratio. Mr. John Pulos of the David Taylor Model Basin recommended a range 
of h/D from 0.05 to 0.10. Also, Timoshenko^ (p.354) has given a value of 
h/D = 0.05 as an estimated division point between thick and thin curved 
bars. Accordingly, the range of values for h/D was selected as 0.05 to 
0 . 10 . 

A calculation revealed that a material with a relatively low modulus 
of elasticity was required in order to cause instability failures in the 
range of h/D greater than 0.05, and, at the same time, remain within the 
upper limit of test pressure of 10,000 psi. Captain Hinners suggested the 
possible use of glass-reinforced plastic, especially in view of the current 
interest in this material because of its favorable strength-weight ratio. 

The David Taylor Model Basin had published a report^® (No. 1413) on 
the hydrostatic pressure tests of a cylindrical shell of a glass fiber 
reinforced epoxy resin fabricated by the Zenith Plastics Co. of Gardena, 
California. The circumferential modulus of this material was reported to 
be 4.8 x 10^ psi, and the plot of pressure vs. strain was linear with but 
a slight departure from linearity near the collapse pressure. This material 
seemed to meet the requirements well. It had disadvantages, however, in 
that it was neither homogeneous nor isotropic, thus violating two basic 
assumptions in the usual approach to the theory of elasticity. 
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It was felt, nevertheless, that these disadvantages would be overcome 
by computing an "effective modulus of elasticity" from strain data, and 
by using models fabricated in such a way as to give maximum dispersion 
of the glass in the resin. Correspondence with the Zenith Plastics Co. 
was initiated to obtain information regarding the procurement of the 
nedessary models. Zenith very generously offered to supply the models at 
no cost. 

At this point specifications were drawn up for the model geometries. 

Using DTMB formula (9) [Formula No. 5, Table 2] , a set of curves was 

drawn for various values of L/D, with buckling pressure as the ordinate 

and h/D as the abscissa. These curves are shown in Fig. 6. The hoop 

stress failure pressure was also plotted on these curves using a nominal 

value of compressive strength of 70,000 psi as obtained from DTMB Report 
38 

1413. These curves were then used to select appropriate values for 
L/D and h/D. It was decided to keep the value of L/D constant at about 
2.5 and vary the value of h/D from about 0.05 to 0.09. It was also 
decided to test a minimum of three models of each geometry in order to 
determine that test results could be reproduced. Table 5 lists the models 
received from Zenith and gives their measured average dimensions. Fig. 7 
is a sketch of the model showing the end closure plugs. The plugs, made 
from medium steel stock, were designed to give simple support of the model 
ends, the condition for the von Mi sea buckling equation. The effective 
length, L, is as shown. The inside diameter was chosen as 3.395" because 
of the availability at the Zenith plant of a mandrel of that diameter. 
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Model Fabrication 



The models were built up using layers of glass fiber reinforced 
epoxy resin tape. The glass fibers ran in one direction, along the 
length of the tape strip, and the tape was applied to the mandrel both 
by circumferential winding under tension (circumferential plies) and 
by laying on strips parallel to the axis of the mandrel (longitudinal 
plies). The process was started and ended with two circumferential plies 
and was continued using alternate layers in a ratio of two circumferential 
plies to one longitudinal ply. This 2:1 ratio of circumferential to long- 
itudinal plies results in the allowable stresses in the longitudinal 
and circumferential directions being in close agreement with the induced 
stresses in these directions. The resulting structure is of course aniso- 
tropic. The thickness was varied by varying the number of layers of tape. 
Five tubes were fabricated, each about 60" in length and each of a different 
thickness (see Table 5). Upon completion of the curing process, most of 
the models were cut from the tube to proper length while still on the 
mandrel at the plant. Some of the models were cut from the remaining 
portions of tubes after delivery, using a band-saw, and grinding or 
machining the ends until smooth and normal to the cylinder axis. By 
weighing and measuring some of the models, the density of the material 
was determined to be about 0.07 pounds/cubic inch. Values of resin 
content for the model material were provided by Zenith, and are given 
in Table 11. 
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FIGURE 7 
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Instrumentation of Models 



Strain data were taken to provide a means for calculating an 
effective modulus of elasticity for the models. As shown in Appendix C, 
values of the strain-pressure ratios, taken in the directions of principal 
stresses, can be used in conjunction with Hooke's Law and thick shell 
theory to solve for and A value for Poisson's ratio must be 

i 

known or assumed. 

In order to determine average values for and strain gages 
were applied to a number of the models, oriented both circumferentially 
and longitudinally. The number of gages used was limited both by 
financial considerations and by the number of leads that could be fed 
through the pressure fittings in the tank. In general, four gages were 
applied to each instrumented model. These gages were located at the mid- 
length of the model to avoid end effects, and were oriented two in a 
circumferential direction and two in a longitudinal direction. The two 
circumferential gages were placed "back to back", inside and outside, as 
were the two longitudinal gages. 

The choice of gage types was governed by: 

(1) the need to have non-pressure-sensitive gages; 

(2) the necessity for measuring strains up to 15,000 microinches/inch ; 

(3) the desire to use a gage length long enough to reduce the effect 
of local variations in strain; and 

(4) the cost of gages. 

On the basis of these considerations, two types of gages were selected. 
For use on the inside of the model, where the gage would not be subject to 
the pressure of the fluid, the SR-4, A-5-1 wire gage was used. This gage 
has a 1/2" gage lenpth, is relatively inexpensive, and is suitable for the 
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measurement of strains up to l£ to 2%. For use on the outside of the 
model, where the gage would be subject to external pressure of the fluid, 
the SR-4, paper-backed conatantan foil gage FAP-50-12 was selected. This 
gage, with a l/2 n gage length, is not as pressure-sensitive as the A-5-1 
gage. It also is suitable for measurements of strain up to 1^- to 2 %, 

Compensating gages of the same type and lot number were used in con- 
junction with the active gages. 

The gages were applied using Duco cement, and the gage manufacturer’s 
instructions were followed. After the cement had dried, all gages were 
connected in the strain indicator circuit and pressed with a pencil eraser 
tip while watching the strain indicator to check for excessive deflections 
of the needle. 
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Teat Procedure 



Strain gages were applied to the inside of the rodel by using a 
"harness"— a sheet of graph paper to which the gages were attached. This 
sheet was trimmed to exact dimensions, rolled into a cylinder, and positioned 
inside the model with glue applied to the gages. Leads of #28 stranded 
wire were attached to the strain gages. 

The ends of the model were prepared by thorough cleaning with acetone. 
PR-1321 Class A sealant * was then applied to the ends, and the end plugs 
inserted, leading the wires from the inside of the model out through the 
1/4" tube *« screwed into the lower end plug. After the sealant set up, 
the model was ready for installation in the tank. Fig. 8 shows a model 
ready for installation and Fig. 5 shows a model in the tank and partially 
lowered into testing position. Fig. 3 shows a model ‘inside the chamber. 

The model was positioned in the open tank by first leading the wires 
from the l/4" end plug tube through the "A" fitting in the center of the 
tank bottom, then connecting the external strain gage leads to the leads 
inside the tank from the "B" fitting. The model was then lowered, allowing 
the l/4" tube to slide through the "A" fitting until the lower end plug 
rested on a steel "spool" support about 2" high located in the bottom of 
the chamber. The lock nut of the "A" fitting was then tightened, taking 
care not to over- tighten. 



* Manufactured by the Products Research Co., 410 Jersey Ave,, Gloucester 
City, New Jersey 

** This tube is a semi -permanent attachment to the end plug, its threaded 
end having been coated with sealant prior to screwing into the tapped 
end plug. 

/ 

*** It was observed that excessive tightening of the lock nut permanently 
deformed the ferrule and tube, making it impossible to slide the tube 
back out of the fitting after the test. 
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With the model in place, the strain gage leads coming from the 1/4" 
tube and from the "B" fitting underneath the pressure vessel were connected 
to the strain indicator or switching unit. The chamber was filled with trans- 
former oil * to a level slightly above the sealing surface. The circular 
plate was inserted, causing the excess oil to be squeezed into the plug 
recess. This precaution was taken to insure that the oil was relatively 
air-free. The screw plug was then positioned and tightened. Tapping the 
screw plug wrench handle with a hammer insured firm seating of the circular 
plate. 

Pressure was raised by operating the hand pump. The valve "G" in the 
line to the 5000 psi gage was shut at 4500 psi. This valve, although rated 

at 3000 psi, performed satisfactorily. 

In some cases, only one strain gage was measured per run, using the 
strain indicator directly. In most of these tests, at least two runs were 
made per gage, and the readings for the two runs were found to vary very 
little. For other tests, the switching unit was employed and a number of 
gage 8 were read per run* 

Upon completion of a test, the screw plug was removed, and the "A" 
fitting lock nut was loosened to permit removal of the model. As the model 
was removed, the l/4" tube slipped out of the "A" fitting, allowing the oil 
in the tank to drain into a receptacle below. 

The two pressure gages used during the tests were calibrated at the 
Material Laboratory before beginning the tests and at DTMB toward the end 
of testing. Both calibrations were made using dead-weight gage 

* VEMC0 Class C transformer oil was used as the pressurizing fluid, making 
it unnecessary to waterproof strain gages exposed to the fluid. 
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testers. No differences were noted in the gage calibrations. The 
5000 psi gage was adjusted after the second calibration to give gage 
readings closer to the actual pressure. The gage calibration curves 
are shown in Figures 9 and 10. 
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EXPERIMENTAL DATA 



A total of twenty- three models was tested, twenty-one of which were 
tested to failure. Seven of the models were instrumented to measure 
strains in order to determine the moduli of elasticity. Three models 
were designed to fail by yielding of the material. This was done in order 
to calculate the maximum compressive stress which the material would 
sustain and to determine the appearance of models after yield failure. 

Since the calculated circumferential stress for the "D” and "E" models 
was close to the nominal strength of the material, this information was 
used to ascertain whether these models did, in fact, fail by instability 
alone, or whether some yielding might have taken place. 

A summary of the models tested, with the respective collapse pressures, 

is presented in Table 6. The recorded strain data are in Appendix A, and 

the pressure-strain plots are in Appendix B. Table 7 summarizes the 

measured strain-pressure ratios as obtained from the plots. The first 
♦ 

models instrumented were tested with the compensating gage located inside the 
model for inside active gages, and with the compensating gage located in the 
pressurized fluid for external active gages. As shown in Appendix C, it is 
theoretically possible to cancel the effect of radial strain on the outside 
active gages by having the compensating gage mounted on a block which is 
subject to the hydrostatic pressure. Furthermore, any effect of pressure on 
the active gage in the fluid would be compensated for , as well as more 
accurate temperature compensation provided. The later models tested, however, 
had all gages — inside and outside — balanced against a single dummy gage in- 
side the model. This latter procedure was followed because it had been 
determined that the error caused by neglecting the Poisson effect of the 
radial strain on the outside gages was negligible, and therefore the 

additional complication to the testing procedure was unjustified. 
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In order to be able to compare and average the strain data taken under 
the two different methods of compensation, the strain-pressure ratios of 
active gages read with the dummy under pressure were corrected to values 
equivalent to those taken with the dummy gage inside the model. This was 
accomplished by simply adding the value of strain-pressure ratio of a gage 
mounted on a dummy block under pressure, balanced against a dummy gage inside 
a model. All the strain-pressure ratios used in computations, therefore, 
have values which would result from active gages being read against a 
dummy located on the inside of the model, not under pressure. 
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ANALYSIS OF EXPERIMENTAL DATA 



Strain Data Analysis 



The moduli of elasticity were calculated from the strain-pressure 
ratios using the formulas presented in Table 8, The theoretical basis 
for these formulas is shown in Appendix C. A summary for the calcula- 
tions for Ef and ^ is given in Tables 9 and 10. No values were cal- 
culated for the B-model because premature failure of the instrumented 
model in this series prevented the taking of sufficient strain data. 

The values of E^ and Ejj used for the B-models are averages of the values 
for the other models. 

Since available instability formulas apply to isotropic materials, 
it was necessary to develop a method for converting the principal mod- 
uli E^ and Ex into an effective modulus E e . This was accomplished by 
first deriving an expression for the shear modulus of an orthotropic 
material. The derivation followed a method used by Timoshenko^ (pp. 
54-57) in deriving the expression for the isotropic case and is presen- 
ted in Appendix D, The shear modulus thus found is 

G = BxEg 

(Ex + E,)(l +*0 

Note that for = E^ this expression reduces to the isotropic case 

0 = — !_ 

2(1 

In the derivation, Poisson’s ratio was assumed to be constant. As noted 
later, this assumption does not lead to an appreciable error. Then, to 
obtain an effective modulus, the orthotropic expression for G was equated 
to that for the isotropic case as follows: 
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Let Ex = kEj* , then G = - . Equating this to 

(1 + y) 

k 

gives 

2(1 +V) 

Ee = 2E < o — - 

1 + k 



The calculated effective values for the modulus are summarized in Table 11. 

In addition, the authors undertook the derivation of an expression 

for the buckling pressure of an orthotropic cylinder. This derivation is 

shown in Appendix E, and the buckling pressure is given in the form of a 

determinant. The calculation of buckling pressure using this determinant 

is very cumbersome . Since the expression was derived shortly before the 

completion of the thesis, no evaluation of the determinant was made with 

respect to the experimental results contained herein. It is believed, 

however, that it may have merit in computing buckling pressures for ortho- 

tropic cylinders. Also, it might be used to verify the assumption of an 

effective modulus for use with the isotropic equations. It is noted that 

the determinant derived for the orthotropic case reduces to the expression 

for the infinite tube, with E^ the only modulus, when L is made to approach 

infinity. This suggests that, in reality, an effective modulus for use in 

the buckling equation should depend on the length of the cylinder. It is 

also noted that the determinant reduces to that given by Timoshenko^® 

(p. 449) when E* is made equal to E* and G is put equal to , 

2(1 

. 2 - 

except that the term - 2- — is added to the last term in row 3, 

h 2 

column C to take care of the end thrust. 

Details of the assumptions employed in the calculation of an effec- 
tive modulus follow. 

(1) As previously stated, Poisson's ratio was assumed to be constant. 
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On the basis of data provided by the David Taylor Model Basin 
and the Zenith Plastics Company, this value was assumed to be 
0.15. The effect of the variation in Poisson's ratio on the 
value of the calculated modulus is small. For example, using 
V - 0.10 instead of 0.15 in the calculation for E^ raises the 
calculated value about 3% for both the A models and the E models. 

(2) The stress field in the model was assumed to be in accordance 
with the Lame thick shell theory, and the strain-pressure 
curves were assumed to be linear up to the point of collapse. 

As a consequence of this assumption the bending stress due to 
the simply supported tube of finite length is neglected. As 
will be shown later, the bending stress prior to buckling is 
a very small quantity. 

(3) The effect of radial strain was neglected in calculating E^ z 

and E Xj . This effect is small as is shown below: 

The radial strain on the outer surface is given by: (see Appendix C) 




Thus, 




It is seen that the greatest error would occur for the largest 



value of m rz , which value would correspond to the thickest 
cylinder, i.e., the cylinder with the least value of Kj. For 



the E-model, ^ = 6.5; K 2 = K x - 1 = 5.5; E^ = 5.5 X 10 6 ; 
E x = 4 x 10^; V - 0.15; then assuming a nominal value for 




Now the exact formula for computing E using outside strain 
gage data is •. t/ . // -t- */V i _ ? 




It was found that the difference between the value of 
calculated by using this equation with = 0.046 and 
that calculated by neglecting m^ entirely was less than 1%. 

It was further noted that the longitudinal compressive strain was 
greater on the outside than on the inside of the models. This difference 
is attributed to the fact that some bending of the shell undoubtedly 
occurs even for the lower pressures. That this difference is not inherent 
in the strain relations can be shown by the following development: 

From Appendix C, 
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*2 “ 




£ 



~pK, 

x 1 “ 2 Ex 




(radial stress = ) 



(radial stress = 0 ) 



Since K 2 = Kj-1, 



€ -€ = 
fc X / - 



E* 




-p _ TS 7° 

~ E* E 



_ -7° AT, , 

Er 2 Ex Et 



from which, 

= v(j=^ - ~j=r) 

Since E^ ^ 5 x 10 6 E^ = 3 x 10 6 and ^=0.15, 
m X2 - b Xi = (0.15) (1/5 - 1/3) = -0.02 
This difference is small in comparison with the differences actually 
measured and is of opposite sign. It is concluded, therefore, that the 
measured difference results from bending of the generators. In order to 
reduce the errors arising from this effect, an average of m x , and m* z 
was used in all calculations which called for either value. 
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Table 7 shows that even for identical models there were considerable 



differences in measured values of strain-pressure ratios. These differences 
may be attributed to two causes; 

(1) The properties of the individual models may vary slightly, even 
though cut from the same tube. 

(2) Even on the same model, the strain-pressure ratios will vary froqa 
point to point around the circumference, as shown in reference (38). 

The most accurate data for use in computing moduli would be that taken 
from a number of strain gages spaced around the circumference of each 
model, using the average values of strain-pressure ratios. Due to limi- 
tations discussed previously, this procedure could not be followed. 

Instead, the averaging process was accomplished by using inside and outside 
strain data, and averaging the moduli thus calculated. Where circumferential 
variations in strain-pressure ratios are caused by the tendency to lobe 
formation, the use of inside and outside data from "back-to-back" gages 
tends to cancel the variation. In the case of the A, D, and E geometries, 
two instrumented models of each were tested in order to provide more 
data for averaging. 

Since the tubes were fabricated using the same techniques and material, 
one would expect the moduli to be about the same for all the models. Such 
is the case, except that the A models were found to have a slightly higher 
modulus. (See Tables 10 and 11). As previously noted, the effective modulus 
for the B models was taken as the average of that of the other models, due 
to insufficient strain data. The small differences in resin content for 
the models (Table 11) are not sufficient to account for the differences 
in -oduli. 
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TABLE 8 



SUMMARY OF THICK SHELL FORMULAS FOR CALCULATING 
MODULI OF ELASTICITY 



Inside gages, dummy block not subject to pressures 



= 



K, (,-**■) 



K, 

Z-trn*, +z/n)p,) 



Outside gages, dummy block not subject to pressure: 

£ _ X 2 . ( nv)(i-2y) Xi ( 

' 7 )4t(t-z') + z'(ynxi+fti*) 2 [>v x z(i-v)+ 



Outside gages, dummy block subject to hydrostatic pressure**: 



(Kz-lXi-v 1 ) 

?n'*z 



r (-*£-- 1 ) 6 -»') 

E-%~ . / 

rn Xa . + v yn?3. 



** Dummy block assumed to be of same material as model; circumferential 
active gages balanced against circumferential dummy and longitudinal 
active gages balanced against longitudinal dummy. 



f 
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TABLE 8 (Continued) 



SUMMARY OF THICK SHELL FORMULAS FOR CALCULATING 
MODULI OF ELASTICITY 




E<p - Circumferential modulus of elasticity 
Ex - Longitudinal modulus of elasticity 



R = 


Mean radius, '/z(R, + R.S) 






■A = 


Wall thickness, (/?*.-£,) 






K t = 


zRl 






Ri-Rf 






Kz = 


Ri+R? 

Rl-Rf- 






m*, - 


Avg. strain-pressure ratio, 
dummy not under pressure. 


inside 


circumferential gages, 


It 

X 

£ 


Avg. strain-pressure ratio, 
dummy not under pressure. 


inside longitudinal gages, 




Avg. strain-pressure ratio, 
dummy not under pressure. 


outside 


circumferential gages 


= 


Avg, strain-pressure ratio, 
dummy not under pressure. 


outside 


longitudinal gages, 


= 


Avg. strain-pressure ratio, 


radial, 


on outer surface. 


II 

5 


Avg. strain-pressure ratio, outside 
dummy under hydrostatic pressure. 


circumferential gages 


II 

v. X 
£ 


Avg. strain-pressure ratio, outside 
dummy under hydrostatic pressure. 


longitudinal gages, 


2V = 


Poisson's ratio. 
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Stress Analysis 



In order to properly evaluate experimental results, it was necessary 
to verify the mode of failure for each model. Although all models were 
designed to fail by instability, the predicted instability failure pressure 
for the D and E models was fairly close to the corresponding compressive 
yield failure pressure predicted using the thick-shell theory and a nom- 
inal value for the yield strength of the material. Also, the thick-shell 
formulas used to calculate the stress do not take into account the bend- 
ing effects. Because of these considerations, the following steps were 
taken in order to make an evaluation for the mode of failure. 

(1) In order to evaluate the effects of bending, a formula was de- 
rived to calculate the longitudinal bending stress for a simply 
supported cylinder. The method of derivation was suggested by 
Mr. John Pulos of DTMB and is presented in Appendix F. The ex- 
pression for the bending stress at midbay is: 






S//vh e £ s/n e 






cos*- e^-hsS/M#*ep, 



\ 



Substituting values for the E model, the theoretical longitud- 
inal bending stress at mid bay was found to be only 75 psi. 

(2) The Lame (thick-shell) stress corresponding to the failure pres- 
sure was calculated for each of the models and compared both 
with a nominal compressive strength *of 70,000 psi as obtained 

33 

from DTMB Report 1413 > and with the lowest calculated value 

obtained from a yield model (63,233 psi). For the E model, the 
highest calculated stress was 59,285 psi (see Table 6). This 
value is 10,715 psi below the value for nominal strength but 
only 3,948 psi below that obtained from the yield model. It 



* Yield strength in circumferential direction, . 



should be noted, however, that application of the Hencky-von 
Miaes failure criterion would cause the differences to be greater. 

(3) Finally, the physical appearance of the collapsed models was 

compared to that of the thinnest model, for which a pure in- 

stability failure was certain, and with that of a short model 
which definitely failed by yielding of the material. Figures 11 
and 12 show collapsed C, D, and E models together with a short 
model which failed by yielding. It is noted that the longer 
models pictured failed along a longitudinal line parallel to 
the model axis by buckling inward, whereas the short model tended 
toward failure throughout most of its wall area. 

On the basis of these considerations it was determined that the A, 

B, and C models failed purely in the instability mode. Since the com- 
pressive stress for the E models was within about 4000 psi of the value 

obtained from a yield model, there still may be reason to question whether 
or not some yielding took place. Also, there was scatter in the points 
for the D and E models (See Fig. 13) which may be a result of local yielding 
caused by stress concentrations due to the heterogeneous nature of the 
material. 
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INDENTATION AND DISCUSSION Df RESin TO 



Resuats of the experiments ere shown graphically in Figure 13. The 

critics! pressures detained fro. tests were converted to the non-dimensional 
form p/g. . These value. plotted againat t(]3 thl0knes3 _ dlameter 

ratio. Calculations for p/E. are shown in TUble 12. 

Also plotted on Figure 13 are curves representing various theoretical 

buckling formulas, using an average value of 1/D for the models tested. 

(Average L/D = 2.53). 

In spite of the anisotropic nature cf the model material and the 

inaccuracies inherent in the computation for an effective modulus, very 

good agreement with theory was obtained emmptfor the thicker models. The 

fact that th, points for the D and K models fall below the curves may indicate 

that the thin shell buckling theory is not adequate for ratios cf h/D in 

this range, is previously noted, however, some yielding may have occurred 
in these models. 

It appears that the von Mises equation is the best one to use in the 

ge f h/D> 0.05 since, in general, it yields conservative values for 

the collapse pressure, assign curves for the buckling cf simply supported 

cylinders were therefore drawn up for the von Hiss. equation. (See Fig. U) . 

Data for these curves were calculated usinv the trm vnon 

mg tne IBM 7090 computer at DTMB. 

These data are reproduced in Appendix G. 

Since curves were not available for determining the number of circum- 
ferential lobes into which a relatively thick cylinder would fall, curves 
were drawn. (See Fig. 15). These curves were calculated using data in 
Appendix 0 by putting the buckling equation in the form 
P/E = P 2 (h/D)3 + F 2 (h/D) 

for each n and L/D and then solving for the desired intersections, 
curves are for v - 0.15. 



These 



For values of n with v - 0.3# in the range h/Di 0.02, see reference 



(33). 

As far as is known by the authors, these curves for the von Mises 
equation and for the number of lobes are the only ones available in the 
range of h/D greater than 0.05. 
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CONCLUSIONS 



Based on the results of these experiments the following conclusions 
are drawn: 

(1) The von Mises equation, equation (6), Ref. 22 is the best in- 
stability equation for cylinders with ratios of h/D greater than 
0.05. It gives conservative values for the buckling pressure up 
to an h/D of about 0.07. Use of presently available simplified 
forms of buckling equations cannot be expected to give conserva- 
tive values for the buckling pressure in this range of h/D. It 
is noted that Vindenburg and Trilling^ had previously reached 
the conclusion that the von Mises formula was probably the best 
one; however, their conclusion was based on an analysis of exper- 
imental and theoretical results for cylinders with ratios of h/D 
less than 0.007 in contrast with the much thicker range of 0.05 
to 0.088 investigated in this thesis. 

(2) For practical cases, the use of stiffening rings, heavy webs, and 
bulkheads would result in smaller ratios of h/D and L/D for a 

given weight of structural material. It is concluded, then, that for 
most practical cases, the thin shell theory is sufficiently 
accurate for predicting instability failures. 

(3) For cylinders built up of fiber glass reinforced plastic layers, 
the methods employed in this thesis are sufficiently accurate 
to predict instability failure pressures. 

(4) For an orthotropic cylinder, the effective modulus, as used in 
the buckling equation, is probably a function of the length 
between supports. 
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RECOMMENDATIONS 



Further Investigations. 

The experimental apparatus designed and assembled for this thesis 
provides a means for the undertaking of a variety of experimental in- 
vestigations in the field of hydrostatically loaded shells. Some investi- 
gations which might be undertaken ares 

(1) Investigate the stresses in short, relatively thick, fiber- 
glass reinforced cylindrical shells. Measure the strain 
distribution along a generator and determine the effects of 
bending. Verify the formula for bending stress which is 
derived in Appendix F. 

(2) Using fiber-glass reinforced plastic models, verify the ortho- 
tropic buckling determinant as derived in Appendix E. In 
conjunction with this, evaluate the determinant to arrive at 
a general expression for buckling pressure similar to the 
von Mises equation for the isotropic case. Using the result, 
develop a method of computing an "effective modulus" for use in 
isotropic buckling formulas. 

(3) Perform experiments on models of unusual form, such as: 

a. Stiffened spheres 

b. Ellipsoids 

c. Sphere-cylinder intersections 



- (>S- 



RECOMMENDATIONS (Continued) 



Improvement of Apparatus. 

Replace valve (G) , (Fig. 4) with one rated at 10,000 psi or higher, 
such as American Instrument Company valve, Cat. No. 44-1505. 

Install a fitting in the line to each gage to damp out the sudden 
drop in pressure when a model breaks, thus protecting the gage pointer. 
There are commercial fittings available for this purpose. 

If a greater number of strain gage leads are required to be run for 
the gages mounted on the outside of the model, it is recommended that a 
C0NAX Corporation thermocouple fitting, Cat. No. TG-20-B8, be installed 
in place of the fitting "A" now installed. The TG-20-B8 will take eight 
#20 wires. 
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APPENDIX A 



RECORD OF STRAIN GAGE DATA 



The following tables consist of the recorded strains corresponding 
to the gage pressures and corrected pressures. The model number, gage type, 
gage factor, date of test, gage location , cycle number, run number, and 
pertinent remarks are also indicated. 

Three types of Baldwin-Lima-Hamilton gages were used for the strain 
measurements: 

Gage Type Nominal Gage Length 

A-5-1 1/2" 

FAP-50-12 1/2" 

FA-100-12 1" 



The type A-5-1 gage was used on the inside of the model only, since 
the wire gage is particularly sensitive to hydrostatic pressure. The 
FAP-50-12, a foil gage, was used both inside and outside the models. The 
FA-100-12 was used on the outside of model A-4 only. All three types 
appeared to give reliable data. No particular difficulty was experienced 
in handling and applying the gages. 

In the row entitled Gage Location . "L" stands for a longitudinally 
mounted gage, "C" for a circumferentially mounted gage, "1" for an inside 
gage, and "2" for an outside gage. 

In some cases the gages were read one at a time requiring a differ- 
ent pressure cycle for each gage. A record was kept of the total number 
of cycles for each model and the number of runs for a particular gage. 
These data were taken in order to ascertain the effects of creep, if any. 
It appeared from the results, however, that creep was not an important 
factor for these short term tests. 
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APPENDIX B 



PLOTS OF PRESSURE VERSUS STRAIN 

This section contains plots of the measured strains. The resulting 
slopes, designated by the letter m with appropriate subscripts are indi- 
cated on the plots. 

It is noted that all plots are very nearly linear until the col- 
lapse pressure is approached. The reversed slope for gage L-2, model 
A-4 indicates the development of a tensile stress as a result of bend- 
ing near the collapse pressure. 
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APPENDIX C 



THEORY FOR CALCULATION OF AND E* FROM STRAIN DATA 



The circumferential stress in a thick-walled cylinder under external 
hydrostatic pressure, p, is given by: 



, . 

* Rl-R? 



\_£er.(4Z)] 



ro rue POINT 



where r is the distance from the center A at which the stress is measured. 
If r is given values of r = R^ and r = R2, the following expressions 
evolve for the circumferential stress on the inner and outer walls of 
the cylinder: 




= 





') 



( 0^ and (J^are positive if compressive.) 

The longitudinal stress is constant over the cross section, and 
is given by: 



r J&- 

< r Ri-R, x 



Letting K, = , 



we have 



(Ref. 42) 



1 > K ' 

T>( K,-l) *= y>Kx 
CTy = r K ‘/z 



(0 



If we assume a material which has three principal values of Young's 
Modulus and a single value of z then from Hooke's Law: 
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Combining Hooke 'a Law and thick shell theory, we have strains in the 

inside and outside surfaces of a thick-walled cylinder subject to external 
hydrostatic pressure; 
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(Note that 0^= p on the outside surface) 
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The stresses in the surfaces of a block of material subject to 
external pressure, p, are given by: 

T? = CT X = CTy - — "p 

From Hooke's Law, the strains on the outside surface of the block 
are given by: 
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The v subscript on a term indicates that it refers to the block rather 
than to the cylinder. 
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Consider the case of strain measurements made with strain gages 
on the inside surface of the cylinder, using a dummy gage not subject to 



pressure. The strains as read are given by* 




ze*) 



Solving the first two of these equations: 
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Where at>*, - e *'/V 

= e *‘/r 



Therefore a model instrumented thus can be subjected to hydrostatic 



an assumed or known value of z' can be substituted to solve for E^. and E^. 
The radial strain does not enter the solution, and the radial stress is 
zero. 

Consider the case of strain measurements made with strain gages 
on the outside surface of the cylinder, using a dummy gage not subject 



pressure and the measured strain-pressure ratios, m and , with 
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to pressure. The strains as read are given by: 
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Solving this system of equations simultaneously gives: 
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Consider the case of strain measurements made with strain gages 
on the outside surface of the cylinder, using dummy gages mounted on a 
block subject to the hydrostatic pressure, the longitudinal active gage 
being balanced against a longitudinally oriented dummy, and the circum- 
ferential active gage being balanced against a circumferentially oriented 
dummy. Assume the dummy block to be of the same material as the cylinder. 

The indicated strain readings ares 
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where the primes indicate that the strains are the Indicated values rather 
than the actual strains. From (3) and (5): 



- / 


- so/ 


/K JL 


z> 


z^AT/ \ 


/ / 


ts 


■ 


- />(- 


£ t 




■ 2£x ; - 


Eft 


Er 


ei 3 


- ~p( 


~LS 


_^/a\ _ 


rD 1 


/ / 


-i; 


x z 


\2-£x 


EL e 


£0 y 


! ( 




' £r 


f 

(z p. 


— -p!~ 


1 


z;K, _ 




£ 


t ( 




r 


T[ 'l 


EV 


2- B * 


E# J 


.£V 


E,. 


r 


- f 


~0<z - 


0 


~p( K '/i - 








€ <&z 


V 

1 


E, 








f 


= / 


r^/2 


-/)_ 


v(Kz. - i ) 


—1 




\ 








Ep 






( 






£ A 

zs 



* 



; V 



PZ 



0 



(/O) 



6^. = O 



- C4- - 



The terms involving radial strain drop out. A comparison of (10) 
with (3) shows that they are identical, except that in (10) the radial 
terms are missing, and K 2 - 1 and K^/2 - 1 appear in place of K 2 and 
Kj/2, respectively. Thus by having the compensating gages arranged as 
in this case, the values of and E^ may be computed from the indicated 
strains 6^ and and the radial strains do not enter the solution. 

Solving (10) s 
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APPENDIX D 



DERIVATION OF SHEAR MODULUS FOR AN ORTHOTROPIC MATERIAL 

Following a method similar to that used by Timoshenko (Ref. 40, 
pp. 54-57) for developing the shear modulus for the isotrlpic case, a 
formula is derived for the shear modulus of an orthotropic material 
by assuming two principal moduli of elasticity and a constant value 
for Poisson’s ratio. 

Consider an element in pure shear. For this condition or = - <r e • rk e ; 
as is readily apparent from Mohr's circle. The element can be represented 
as shown below: 







(The deformed element is represented by dotted lines.) The deformed 
angle at A is-^r-Y'j that at B is -^r + Y , where Y Is the shearing strain. 
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arrived at a similar expression using an energy analysis. 
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APPENDIX E 



DERIVATION OF EQUATION FOR THE CRITICAL PRESSURE 
OF A SIMPLY SUPPORTED, ORTHOTROPIC, CIRCULAR 
CYLINDER SUBJECTED TO HYDROSTATIC PRESSURE 

Stress-strain relations ♦ 

For isotropic shells, the relation between the stress resultants 
and couples and the corresponding strains is: 



where the subscript '©"refers to strains in the middle surface. 

If is is assumed that the cylinder is orthotropic with principal 
moduli Ex and E© and with a constant value of Poisson's ratio, these 
equations may be written: 




/V\* = ~D(X x + >^X e ) 





N KQ — G-K 



/V\ xe = DC/-*') Xxe 







n k& — G- yz 




- El- 



The strains and changes of curvature in terms of the displacements 
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The stress resultants and couples can then be written in terms of 
the displacements: 
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The Equilibrium Equations, 

The equilibrium equations for a circular cylinder under radial 
load only are: 
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See Timoshenko, Ref. 30, p.447j note that the term fa.\l T7e x / 

has been neglected, as was done by von Mises^. J 

% 

These equations may be modified to include the effect of axial 
load as follows. This procedure was used by von Mises. 

The total axial force on the ends of the cylinder is equal to 7T p 
This force is transmitted to the ends of the cylinder of area zrrcch . 
Hence the pressure in the cylinder walls is 



z _ rraT P _ P 
P ~ Z 7 T 6 *h ~~ Zh 



Now consider an element of the shell in the longitudinal plane. 
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After deformation, the curvature of the generatrix (neglecting 



higher order terms) is 

, z , 
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T 7 1 



The force per unit length, p'h, thus has a component in the a-direction 



of 



p h ~ Jx l 



Since the radius, a, had been divided out of the original equilibrium 
equation, this term is written as 
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In order to include the effect of the axial load, this term is added 
to the term in parenthesis in equation (3) to getj 

y/iK _l y^Ae / y^r 

+ A d e*- + No - P \ + jj-f. +■ 2. ) — O 

Boundary conditions. 

The transverse and longitudinal cross seotions of the cylinder assume 
the following shapes after buckling j 
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At x = 0, the displacements may be written: 



U = A s/n n & 
s\r — B cos n 0 

■us — C sin n e 



*\T /s 30° out of ph^se. cutth U «W. 



At x = 0, u = 0 and v and w are maximum; hence, the displacements for 
any value of x become: 



u — A s/n n o s/m jr 

a „ ns 

/i r — d cos ns cos £ 

^ rrx 

aj/~ — C s/m h O cos —ji 
ira. 

Let \ ~ -JC t then 

U = A s/n n e s/n /ii. * 
/\ r = d cos ne cos */a. K 
Ajs = C s/m no Cos 



Determinant for the Critical Pressure, 

By substituting the boundary conditions into the relations for the 
stress resultants and couples, and substituting the resulting expressions 
into the equilibrium equations, a system of three equations in A, B, and 
C is determined. In order that an instability condition exist, the constants 
A, B, and C, which determine the magnitudes of the displacements, must 
have values other than zero. This condition requires that the determinant 
of the coefficients of A, B, and C be zero. The determinant for this case 
is given on the next page. Evaluation of the determinant for a given 
geometry and material yields the critical buckling pressure. 
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APPENDIX F 



DERIVATION OF AN EXPRESSION FOR THE LONGITUDINAL BENDING STRESS 

IN A 

SIMPLY SUPPORTED CIRCULAR CYLINDER SUBJECTED TO EXTERNAL HYDROSTATIC 

PRESSURE 



Salerno and Pulos^3 have derived expressions for the stresses in 
a circular cylindrical shell supported by uniformly spaced circular 
rings of constant cross section. By using the displacement function 
obtained in their paper with the appropriate boundary conditions for 
simple support, an expression for the longitudinal bending stress at 
midbay is derived. 

The longitudinal bending stress is given by: 

_ Eh 

^ ~ 2.(1- ■**■) <L X* 



The displacement function, xv-, as derived by Salerno and Pulos for 
deflection in the radial direction, is 



sus = B Cosh (A, *) + F cosh (As x) — 
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where B and F are constants depending on the boundary conditions, A , 7*= A 
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Taking the x-axis along the axis of the cylinder, with the origin at 
the mid-length, the boundary conditions for simply-supported ends require 



that 



from which, 
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o - B A ^ cosh A , Vz +FAS COSH Aj Vz. 



Solving simultaneously, 
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Krenzke and Short ^ have presented a graphical method for calculating 
the stresses obtained by the analysis of Salerno and Pulos. In this 
reference, different parameters were employed as a simplification. 

These parameters, defined below, are used for the remainder of this 
derivation, 
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Using the Krenzke and Short parameters, 
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Taking the second derivative of this expression, simplifying, and sub- 
stituting into the expression for the longitudinal bending stress gives: 
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APPENDIX G 



SOLUTION OF VON MISES 1 EQUATION (6) 

This appendix contains solutions of the von Mises equation (6)^2 
as computed by the IBM 7090 computer at DTMB. Values of p/E were com- 
puted for a range of L/D from 0.50 to 3.00 in increments of 0.25, and 
for a range of h/D from 0.01 to 0.10 in increments of 0.01. The solution 
for a particular h/D was started with a value of n = 2. The value of n 
was then incremented by one until a minimum value for p/E was obtained. 
Columns 3, 4, and 5 give the auxiliary quantities and for 

use in the von Mises equation. 
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